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1.0.  INTRODUCTION 

Since  1980,  the  Applied  Research  Function  has  been  investigating  advanced 
suspension  system  concepts  for  the  Track  and  Suspension  Sub-Function.  These 
concepts  all  presume  on-the-move  adjustability  in  the  springing  and/or  the 
damping  characteristics  of  the  suspension  system.  In  order  to  design  this 
type  of  suspension,  it  was  recognized  from  the  beginning  that: 

o  a  system  would  be  required  for  sensing  some  of  the  vehicle  states, 
and  perhaps  for  sensing  and  interpreting  the  terrain; 

o  that  control  logic  would  be  needed  for  determining  the  optimum 
suspension  adjustments  as  a  function  of  terrain  and  vehicle 
parameters ; 

o  that  an  on-board  microprocessor  would  be  needed  to  carry  out  the 
required  computations  and  for  actuating  the  controllable  suspension 
elements* 

Thus,  the  overall  approach  was  to  identify  in  the  technical  literature  the 
current  state-of-the-art  approaches  to  optimal  suspension  design  and  control 
strategies;  evaluate  these  strategies  for  potential  application  to  combat 
vehicles;  develop  new  strategies  or  modify  existing  ones  to  best  suit  the 
program  objectives;  and  simulate  selected  strategies  (either  found  in  the 
literature  or  developed  by  the  investigators)  in  a  validated  vehicle  ride 
quality  computer  model. 

Section  5.0.  of  this  report  documents  in  historical  order  the  investigations 
conducted  and  the  improvements  in  ride  quality  achieved  with  advanced 
suspension  systems  designs  and  control  strategies. 


2.0.  OBJECTIVE 

The  primary  objectives  of  this  advanced  suspension  technology  investigation 
are: 

o  to  identify  and  simulate  vehicle  suspension  system  control  strategies 
which  will  extend  combat  vehicle  operational  capabilities  to  higher 
speeds  and  rougher  terrains  by  increasing  the  overall  comfort  of  the 
crew,  the  tracking  capability  of  the  gunner,  and  the  vibrational 
environment  of  the  fire  control  system; 

o  to  specify  and  identify  potential  controllable  suspension  system  com¬ 
ponents,  sensors,  and  microprocessor  hardware  required  to  implement 
and  field  test  a  suspension  system  control  strategy. 


3.0.  CONCLUSIONS 

Terrain  elevation  sensing  on  board  a  moving  vehicle  could  be  accomplished  bv 
a  sonar  reflective  sensor  or  a  laser  used  in  combination  with  wheel-mounted 
accelerometers.  (Section  5.1.  and  Appendix  A). 
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Digital  filtering  and  other  digital  signal  processing  techniques  using  the 
1 6-bit  arithmetic  available  with  current  microprocessors  can  be  used  to 
identify  the  overall  roughness  and  the  wavelength  composition  of  a  terrain 
from  a  series  of  terrain  elevation  values.  (Section  5.2.  and  5.'^.). 

The  optimum  damping  level  for  a  given  vehicle  on  a  particular  terrain 
segment  can  be  well-estimated  by  a  combination  of  linear  regression  and  the 
computation  of  a  frequency  normalized  terrain  roughness  calculation  tech- 
nique.  (Section  5«4«). 

Substantial  improvement  in  comfort-limited  vehicle  speed  is  attainable  with 
a  relatively  simple  scheme  which  adjusts  damping  level  as  a  function  of 
time-averaged  vehicle  pitch  rate.  (Section  5.5.). 

Dramatic  improvements  in  vehicle  speed  and  the  elimination  of  resonant 
response  can  be  obtained  by  using  passive  on-off  dampers  with  the 
appropriate  control  logic.  (Section  5.6.). 


U . 0 .  RECOMMENDATIONS 

Extension  of  the  investigation  on  advanced  suspension  technology  concept  may 
include  the  following: 

Further  investigation  on  the  digital  filtering  techniques  proposed  in  this 
report  for  determining  optimum  adaptive  damping  levels  for  vehicle  suspen¬ 
sion  system.  The  developed  suspension  control  law  may  be  tested  using  a 
laboratory  adaptive  suspension  concept  model  or  through  simulations. 

Development  of  adaptive  suspension  control  law  using  direct  parameter  opti¬ 
mization  technique.  The  optimization  objective  may  consist  directly  of 
minimizing  the  average  power  absorbed  at  a  given  location  on  a  vehicle.  The 
resulting  control  algorithm  may  Involve  on-off  (bang-bang)  control  logic  and 
continuous  control  logic* 

Development  of  adaptive  suspension  control  logic  based  on  overall  vehicle 
suspension  performance.  A  model  reference  adaptive  control  suspension 
scheme  may  be  suggested. 

Development  of  optimal  suspension  control  system  using  optimal  and  preview 
control  techniques.  The  techniques  consist  of  sensing  or  predicting  future 
inputs  (e.g.,  the  terrain  being  approached)  and  employing  the  Information  to 
control  and  improve  the  vehicle  ride. 

Extensive  digital  computer  .simulation  for  verifying  potential  ride  improve¬ 
ment  that  can  be  achieved  in  military  vehicles  using  the  proposed  adaptive 
suspension  concepts.  Realistic  vehicle  data  and  test  terrain  profiles  may 
be  used  for  conclusive  results. 

Development  of  microcomputer  controller  system  for  real-time  implementation 
of  proposed  suspension  control  logic.  A  microcomputer  in-the-loop  real-time 
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simulator  consisting  of  a  hybrid  computer  military  vehicle  model  interfaced 
to  the  microcomputer  damper  controller  may  be  carried  out  to  experimentally 
verify  the  potential  ride  improvement. 

Construction  of  laboratory  microcomputer  controlled  suspension  concept 
models  for  complete  hardware  in-the-loop  model  simulation.  The  laboratory 
suspension  models  may  be  used  for  demonstrating  and  testing  various  adaptive 
suspension  concepts.  The  experience  acquired  in  one  laboratory  experiment 
will  be  invaluable  for  future  implementation  effort  as  well  as  for 
recognizing  potential  practical  problems. 

Study  and  investigation  on  current  state-of-the-art  sensor  methodology  for 
on-board  measurements  of  vehicle  states  and  terrain  profile. 

Study  and  investigation  on  current  state-of-the-art  suspension  elements  such 
as  heavy-duty  low-power-input  variable  dampers. 

Development  of  in-house  research  capability  for  implementing  and  evaluating 
advanced  adaptive  suspension  concepts.  Close  interaction  with  academic 
institutions  and  industries  may  be  suggested  for  a  quick  buildup  of  such 
capability. 


5.0.  DISCUSSION 

5.1 .  Preliminary  Study 

In  September  1980,  Francis  Hoogterp  issued  a  report  "Preliminary  Study  for 
Adaptive  Suspension  System  Program."  It  is  Included  as  Appendix  A  to  this 
report.  The  first  technical  area  it  addressed  was  the  need  to  sense  terrain 
elevations  at  1-foot  intervals  at  speeds  up  to  50  mph  (73  ft/sec).  It  was 
concluded  that  the  distances  Involved  are  too  short  for  radar,  but  a  sonar 
reflective  sensor  made  by  Polaroid  might  be  modified  to  meet  this  require¬ 
ment.  A  pressure  transducer  located  within  a  hydropneumatic  type  suspension 
Unit  or  an  accelerometer  mounted  on  a  roadarm  near  the  roadwheel  could  a1 so 
be  used  to  sense  terrain  elevations.  For  all  of  these  three  schemes,  a 
hull-mounted  accelerometer  would  also  be  needed  to  remove  the  effect  of 
hull  vertical  motion. 

Assuming  the  terrain  elevation  can  be  sensed  at  1-foot  intervals,  the  next 
step  is  to  compute  the  RMS  and  perhaps  frequency  content  of  the  terrain. 

Mr.  Hoogterp  suggests  that  one  high  pass  and  three  low  pass  digital  filters 
implemented  on  a  1 6-bit  microprocessor  could  be  used  to  meet  this  require¬ 
ment  . 

The  Environmental  Research  Institute  of  Michigan  (ERIM)  is  developing  a 
laser-based  terrain  previewing  system  for  the  Department  of  the  ^rmy 
Advanced  Research  Projects  Agency.  It  will  have  a  40°  field  of  view  in  azi¬ 
muth  and  elevation,  a  ranging  accuracy  of  i3mm  RMS,  operate  on  7mw  of  laser 
power,  have  a  30-foot  look-ahead  capability,  and  fit  in  a  4-lnch  cube.  It 
is  intended  for  use  on  a  slow-moving  hexapod  "walking  machine".  ERIM  seems 
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to  have  the  expertise  to  modify  this  system  for  vehicle  ride  purposes  and  to 
process  the  acquired  data  signals  in  order  to  characterize  the  terrain.  The 
cost  for  a  prototype  unit  would  most  likely  exceed  $100,000. 

5.2.  Terrain  Characterization  Methods  for  Microcomputers 

5.2.1.  Introduction 

Since  a  knowledge  of  the  terrain  would  be  an  Important  input  to  a 
controllable  suspension  system,  digital  computer  experiments  were  conducted 
to  determine  the  potential  of  digital  filters  for  identifying  the  overall 
roughness  and  the  frequency  content  of  terrains. 

This  effort  started  with  the  assumption  that  the  terrain  elevation  could  be 
sensed  at  1-foot  intervals  and  converted  to  a  digital  value  and  made 
available  to  the  computing  algorithm.  The  goal  was  then  to  find  mathemati¬ 
cal  methods,  suitable  for  microprocessor  implementation,  which  can  detect 
both  the  roughness  and  frequency  content  (wavelengths)  of  the  terrain. 

5.2.2.  Digital  Filtering  Techniques 

The  main  advantages  of  digital  filtering  over  the  more  standard  mathematical 
methods  (such  as  Fourier  transforms)  are: 

o  much  smaller  computer  memory  requirements; 

o  much  faster  computation  time  (to  permit  real-time  performance); 

o  updating  of  terrain  parameters  on  a  continuous  (on  the  order  of  once 
per  foot)  basis  rather  than  on  an  occasional  (once  per  250-500  foot) 
basis. 


The  disadvantages  are  some  loss  of  accuracy  in  the  overall  RMS  and  much  1  ess 
detail  in  distinguishing  component  frequencies.  However,  for  the  purpose  of 
this  project,  a  5K  RMS  error  would  be  acceptable  and  it  is  only  necessary  to 
identify  the  one  or  two  regions  of  the  frequency  spectrum  which  typically 
account  for  most  of  the  terrain  roughness. 

Digital  filters  were  designed  to  compute  the  overall  RMS  of  a  terrain  eleva¬ 
tion  input  signal  and  to  compute  RMS  values  of  the  terrain  from  four  dif¬ 
ferent  frequency  regions.  The  frequency  regions  were  defined  by  the  period 
of  the  waveform  being  2  to  10  feet,  20  to  MO  feet,  and  40  to  80  feet.  Waves 
shorter  than  two  feet  cannot  be  detected  with  a  1-foot  sampling  Interval, 
while  waves  longer  than  80  feet  do  not  appreciably  influence  ride  quality 
for  the  range  of  vehicle  speeds  being  considered. 

The  first  filter  design  considered  was  a  system  of  Butterworth  type  high 
pass  and  band  pass  filters  of  fourth  order  implemented  in  60-bit  floating 
point  arithmetic.  The  cut-off  frequency  of  the  high  pass  filter  was  set  at 
1/80  Hz.  in  order  to  filter  out  wavelengths  greater  than  80  feet.  The  root 
mean  square  (RMS)  value  of  the  high  pass  filtered  signal,  Yi,  was  computed 
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where  N=512.  This  was  considered  to  be  the  overall  RMS  of  the  signal  as 
shown  in  Table  5-1. 

The  original  signal  was  also  filtered  using  the  four  pole  band  pass  filter 
with  the  cut-off  frequencies  set  for  80-  and  40-foot  waves.  The  RMS  of  this 
filtered  signal  is  designated  in  Table  5-1  as  the  "80-40"  foot  RMS  of  the 
signal.  By  setting  the  desired  cut-off  frequencies  in  the  band  pass  filter, 
the  same  method  was  used  to  compute  the  "40-20"  foot  RMS  and  the  "20-10" 
foot  RMS  of  the  input  signal.  Finally,  the  high  pass  filter  with  the  cut¬ 
off  frequency  set  at  1/10  Hz.  was  used  to  compute  the  "10-0"  foot  RMS. 

This  system  of  digital  filters  was  tested  to  determine  its  ability  to  iden¬ 
tify  the  frequency  of  a  simple  sine  wave  with  a  il 2-inch  amplitude.  Several 
frequencies  were  input  (one  at  a  time)  near  the  center  and  at  the  edge  of 
the  frequency  bands.  The  results  are  shown  in  Table  5-1. 

Table  5-1.  Four-Pole  Butterworth  Filter  Results  (60-Bit  Arithmetic) 


Input 

Wavelength 

Overall 

RMS 

80-40 

Computed 

40-20 

RMS 

20-10 

10-.0 

4 

8.485 

5.837 

.003 

.008 

8.485 

8 

8.485 

1.996 

.016 

.918 

7.936 

10 

8.496 

.205 

.051 

6.002 

6.002 

14 

8.477 

.554 

.376 

8.480 

2.007 

20 

8.466 

.630 

6.547 

5.386 

.393 

28 

8.478 

.401 

8.521 

.463 

.122 

40 

8.421 

5.964 

5.964 

.712 

.029 

56 

8.277 

8.433 

.454 

.012 

.008 

80 

5.985 

6.019 

.524 

.032 

.0026 

119 

1.672 

.299 

.009 

0 

0 

The  theoretical  RMS  of  each  of  the  signals  is  8.485  inch.  Note  that  for  an 
80-foot  period,  the  overall  RMS  is  only  5.985.  This  is  the  result  of  the 
80-foot  cut-off  frequency  in  the  high  pass  digital  filter,  which  is  expected 
to  reduce  the  amplitude  of  an  input  wave  at  the  cut-off  frequency  by  3db 
(l.e.,  multiply  it  by  ^/2).  The  amplitude  reduction  for  the  119-foot  wave 
is  much  greater  as  shown  by  the  total  RMS  value  of  only  1.672.  Whenever  the 
actual  period  is  in  the  center  of  one  of  the  digital  filter  pass  bands,  the 
signal  frequency  is  very  strongly  identified  with  that  band  alone.  For 
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instance,  for  a  28-foot  input  wavelength,  Table  5-1  shows  an  RMS  value  of 
8.521  in  the  "40-20"  foot  region.  When  it  is  at  the  boundary  of  two  adja¬ 
cent  pass  bands  0.707  of  the  total  signal  RMS  is  assigned  to  each  pass  band. 

If  it  is  desired  to  always  characterize  the  terrain  by  a  single  frequency 
waveform,  then  a  scheme  could  be  readily  devised  to  estimate  the  frequency 
and  amplitude  of  the  wave  from  the  values  produced  by  this  system  of  filters 
(for  example,  as  the  maximum  value  of  a  parabola  computed  to  fit  the  four 
RMS  vs.  frequency  values).  The  "80-40"  foot  band  pass  filter  seems  to  have 
some  problems  with  wavelengths  of  less  than  10  feet,  and  therefore,  requires 
more  investigation.  The  sampling  interval  of  one  foot  may  be  responsible, 
and  it  probably  also  accounts  for  the  scatter  in  the  rest  of  the  table 
entries . 

The  previous  digital  filtering  results  were  achieved  on  a  60-bit  computer. 
Since  digital  filtering  is  known  to  be  sensitive  to  the  number  of  digits 
used  in  the  computations,  these  filters  were  also  tested  on  a  computer 
operating  with  32-bit  floating  point  arithmetic.  The  high  pass  filtering 
was  done  with  a  four-pole  Butterworth  high  pass  digital  filter.  A  four-pole 
Butterwprth  band  pass  filter  was  tried  for  the  required  band  pass  filtering, 
but  it  did  not  function  properly  because  of  the  32-bit  word  size.  This 
problem  was  overcome  by  applying  a  two-pole  Butterworth  band  pass  filter  to 
the  signal  two  times  in  succession. 

As  was  done  with  the  60-btt  filters,  a  number  of  computer  runs  were  made 
with  artifically  constructed  input  signals.  This  was  done  so  that  the  RMS 
values  produced' by  the  program  in  each  of  the  four  frequency  intervals  could 
be  compared  with  the  theoretical  values  of  the  input  signal. 

The  results  as  displayed  in  Table  5-2  show  that  when  only  a  single  frequency 
wave  is  present  and  is  at  the  center  of  one  of  the  four  frequency  intervals, 
it  usually  results  in  the  RMS  value  for  the  proper  frequency  interval  being 
about  nine  times  as  high  as  for  any  other  of  the  frequency  intervals.  When 
the  frequency  of  the  input  wave  is  at  the  boundary  of  two  frequency  inter¬ 
vals,  equal  RMS  values  result  in  these  two  intervals,  and  the  RMS  in  the  two 
remaining  intervals  is  lower  by  more  than  a  factor  of  ten.  For  wavelengths 
of  less  than  80  feet,  the  overall  RMS  of  the  terrain  is  related  approximate- 
ly  to  the  four  frequency  filtered  BMS  values  by: 

RMS2  tot  =  Bo-40  +  RMS2  no-20  +  RMS2  20-10  +  ^MS^  10-0 

This  can  be  seen  in  both  Tables  5-1  and  5-2. 

For  the  four  frequency  intervals  studied  in  Tables  5-1  and  5-2,  let  the 
signal-to-noise-ratio  be  defined  as  the  result  of  dividing  the  calculated 
RMS  value  which  is  closest  to  the  input  wavelengths  by  the  largest  calcu¬ 
lated  RMS  value  among  any  of  the  other  wavelengths.  The  signal  to  noise 
ratios  calculated  from  the  values  displayed  in  Tables  5-1  and  5-2  are  shown 
in  Table  5-3. 

The  overall  results  of  this  section  show  that  this  type  of  filtering  scheme 
can  identify  the  various  wavelengths  present  in  a  terrain,  even  when  a 
32-bit  computer  is  used. 


Special  provisions  were  made  in  the  32-bit  general  purpose  computer  program 
to  simulate  the  effect  of  1 6-bit  integer  microprocessor  arithmetic.  With 
this  scheme,  the  double  application  of  the  four-pole  Butterworth  band  pass 
filters  showed  a  tendency  toward  Instability  and  erratic  results.  However, 
the  two-pole  filters  applied  four  times  to  a  single  frequency  sine  wave  pro¬ 
duced  the  results  shown  in  Table  5-4.  The  RMS  value  at  the  middle  of  the 
frequency  interval  is  within  2%  of  its  theoretical  value.  However,  the  edge 
values  are  3115  lower  than  for  the  four-pole  filter  applied  twice. 
Theoretically,  the  RMS  at  the  edge  should  be  /?/2  times  the  center  RMS. 

This  enables  the  total  RMS  of  the  signal  to  be  computed  from  the  RMS  in  two 
neighboring  intervals.  Thus,  the  edge  values  in  Table  5-4  (for  the  20-  and 
40-foot  wave)  are  about  30^  lower  than  ideal.  However,  the  filter  coeffi¬ 
cients  could  be  adjusted  to  rectify  most  of  this  problem. 

All  the  computations  described  above  have  been  carried  out  on  fixed  sets  of 
numbers  of  length  512.  As  a  vehicle  is  traversing  terrain,  any  such  array 
of  terrain  elevations  would  need  to  be  updated  at  every  sample  interval  by 
adding  the  most  recent  elevation,  deleting  the  oldest,  and  recomputing  the 
RMS  values.  A  more  practical  approach  is  to  submit  each  new  terrain  eleva¬ 
tion  to  the  digital  filters  in  order  to  produce  four  filtered  terrain 
signals.  For  each  of  these  signals,  an  N-point  running  average  mean  square 
(MS)  value  can  then  be  computed  using  the  formula: 

MSt  =  +  (N-1)MSt>J/N 

N  should  be  chosen  consistent  with  the  frequency  of  the  filter  which  pro¬ 
duced  Y  from  the  terrain  elevation.  This  method  produced  instantaneous  RMS 
values  corresponding  to  a  recently  traversed  section  of  terrain  whose  length 
is  determined  by  N.  This  method  was  tested  in  32-bit  arithmetic  and  pro¬ 
duced  accurate  and  non-oscillatory  results  as  long  as  N  was  longer  than  the 
longest  wavelengths  in  the  filtered  terrain  elevation  signal. 

5.2.3.  Running  Average/Predominant  Period  Technique 

The  data  displayed  in  Table  5-5  are  the  result  of  estimating  the  single  pre¬ 
dominant  frequency  in  a  data  signal  by  a  means  other  than  digital  filtering. 
Running  average  values  are  computed  continuously  over  the  most  recent 
50-foot  interval  for  the  RMS  values  of  the  terrain  elevation  (RMSE)  and  the 
terrain  slope  (RMSS).  .An  estimate  for  the  length  or  period,  P,  of  the  most 
predominant  waveform  present  in  the  terrain  is  then  computed  by: 

P  s  2'«’  RMSE /RMSS 

The  development  of  this  formula  is  documented  in  Reference  1.  The  range  of 
values  obtained  in  Table  5-5  is  due  to  the  50-foot  running  average  effect, 
and  perhaps  a  longer  average  should  be  used.  The  mathematical  development 
of  this  method  is  documented  in  Reference  1.  The  running  average  com¬ 
putation  can  be  efficiently  carried  out  by  the  following  scheme: 

(1)  Let  Xi  be  the  ith  data  value. 
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Table  5-2.  Two-Pole  Butterworth  Filter  Applied  Twice  (32-Bit  Arithmetic) 


Input 

Wavelenf?th 

80-40 

Computed  RMS 
40-20 

20-10 

10-0 

Total 

RMS 

•3.5 

.001 

.000 

.008 

8.477 

8.485 

10 

.005 

.123 

6.103 

6.002 

8.496 

14 

.024 

.837 

8.480 

2.007 

8.477 

20 

.133 

6.036 

5.994 

.479 

8.476 

28 

.853 

8.521 

.997 

.122 

8.478 

40 

6.038 

6.033 

.143 

.029 

8.421 

56 

8.428 

.986 

.029 

.008 

8.277 

119 

.670 

.021 

.001 

.000 

1.672 

Table  5-3. 

Digital  Filtering  Signal-to-Noise  Ratios 

Input 

Wavelength 

60-Bit  Arithmetic 
Signal-to-Noise  Ratio 

32-Bit  Arithmetic 
Signal-to-Noise  Ratio 

10 

14 

20 

21 

28 

40 

56 

80 


41 .41 
4.23 

«  ' 

13.46 
18.40 
11 .85 
18.57 
11.49 


69.59 

4.23 

17.76 

« 

8.55 

59.69 

8.55 

» 


•Not  Available 
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I  Table  5-4.  Two-Pole  Butterworth  Band  Pass  Filter  Applied  Four  Times, 

1 6-Bit  Integer  Arithmetic 


Period _  40-20  RMS 


3.5 

.000 

10 

.240 

a 

14 

1.049 

18 

2.962 

20 

4.347 

22  ■ 

5.858 

24 

7.255 

28 

8.383* 

32 

7.464 

36 

5.676 

1 

40 

3.967 

42 

3.283 

56 

.872 

80 

.076 

119 

.000 

*Theoretical  =  8.4385 

_ _ _ _ _ : _ _ _ — - 

TABLE  5-5.  Frequency  Estimation  By  Running  Average  Method 


Actual  Estimated 

Period  ^ _  Period 


8 

8. 1-8.4 

10 

10.0-10.4 

14 

13.8-14.5 

21 

20.4-21.7 

28 

26.8-29.4 

40 

39.2-40.8 

56 

51-61 

80 

72-89 

119 

101-140 
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(2)  Compute  A  =  2,  ^ 

50  S  *1 

i=1 

(3)  For  X51,  compute  A  s  JL  C50A  -  xi  +  X5-1I,  and  store  the  value  of 

50 


X51  in  x-i. 

(4)  For  X52  thru  xioo>  ^^®  ^  range  from  52  to  100  and  com- 

pute:  ,  ^  / 

A  =  J.  t50A  -  Xi_50  +  *11 »  store  the  value  of  Xj^  in  Xj^_5o. 

To 


(5)  For  X101  thru  xi5o>  continue  the  process  as  in  steps  3  and  4, 
but  use  x^oi  thru  X')50  in  place  of  X51  thru  X'^qo* 

These  preliminary  studies  show  that  properly  devised  digital  filtering  and 
other  signal  processing  techniques  can  be  used  to  assess  the  overall  rough¬ 
ness  of  a  terrain  and  its  wavelength  composition.  They  indicate  that  a 
16-bit  digital  microprocessor  should  be  sufficiently  accurate  to  obtain  the 
desired  results.  Actual  implementation  on  microprocessor  hardware  will  be 
required  to  determine  if  the  computations  can  be  done  fast  enough  for  on¬ 
board  vehicle  operations.  This  discussion  in  Section  5.3.  addresses  this 
question. 

5.3.  Design  and  Microprocessor  Implementation  of  Digital  Filters 

Dr.  K.  C.  Cheok  and  Dr.  N.  K.  Loh  of  Oakland  University  have  completed  a 
comprehensive  and  technically  in-depth  study  entitled ,  "Design  and 
Microprocessor  Implementation  of  Frequency  Selective  Digital  Filters  for 
Application  to  Terrain  Roughness  Identification,"  Reference  2.  The  theoreti¬ 
cal  and  mathematical  basis  for  analog  and  digital  signal  filtering  was 
thoroughly  reviewed  and  carefully  applied  to  the  terrain  roughness  problem. 
Technical  problems  such  as  frequency  warping,  distortion,  and  aliasing  can 
occur  because  of  the  discrete  sampling  of  the  signal  and  the  number  of 
digits  (or  bits)  available  in  mloroprooessor  arithmetic.  Methods  to  avoid 
or  overcome  these  problems  were  incorporated  into  the  filter  design  process, 
taking  into  account  the  characteristics  of  the  particular  microprocessor 
to  be  used— the  Motorola  MC6802.  Successful  Implementation  on  the  micropro¬ 
cessor  required  development  of  fast  and  accurate  arithmetic  subroutines,  the 
scaling  of  the  recursive  digital  filter  equations,  the  special  handling  of 
arithmetic  saturation,  and  careful  attention  to  memory  management. 

A  laboratory  set  up  was  created  at  Oakland  University  in  order  to  observe 
the  actual  performance  of  the  total  digital  filtering  system.  The  experi¬ 
mental  results  generally  agree  well  with  the  theoretical  predictions,  though 
there  are  minor  discrepancies  which  can  be  remedied  through  further  fine 
adjustments  in  the  set-up.  The  system  demonstrates  sufficiently  fast  com¬ 
putational  speed  and  adequate  accuracy  for  use  in  the  identification  of 
terrain  frequencies. 
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5.4.  Optimum  Constant  Damping  Level  as  a  Function  of  Terrain 

It  has  long  been  observed  from  both  vehicle  field  testing  and  computer  simu¬ 
lation  that  a  higher  level  of  damping  can  improve  vehicle  ride  on  rougher 
terrains  while  giving  a  poorer  ride  on  smoother  terrains,  and  vice-versa. 
Thus,  vehicle  damping  design  is  usually  a  compromise  aimed  at  moderate  or 
average  terrain  roughness  conditions.  One  of  our  major  goals  was  to  find  a 
method  for  predicting  optimum  damping  as  a  function  of  terrain  roughness 
and/or  frequency  content  in  a  way  which  could  be  Implemented  in  real-time 
using  an  on-board  microcomputer. 

One  method  of  accomplishing  this  goal  is  documented  in  this  section.  It 
involves  simulating  the  LK10934  (a  20-ton  tracked  combat  vehicle  concept) 
over  nine  different  terrains  using  the  TACOM-validated  VRIDE  computer  model, 
Reference  3.  On  each  terrain  the  level  of  damping  (or  damping  factor),  was 
determined  which  allowed  the  vehicle  to  achieve  the  highest  speed  at  the 
ride-comfort  limiting  level  of  6  watts  driver  vertical  absorbed  power.  This 
method  of  measuring  ride  comfort  is  in  Army-wide  use  and  is  documented  in 
Reference  4.  Damping  was  varied  by  multiplying  the  entire  force  vs.  velo¬ 
city  curve  by  a  constant  factor.  A  given  level  of  damping  and  a  constant 
vehicle  speed  were  maintained  for  each  run  over  each  course. 

Terrain  RMS  roughness  varied  from  0.67  to  3.88  inches,  vehicle  speeds  ranged 
from  15  to  85  mph,  and  the  optimum  damping  level  varied  from  0.4  to  1.4 
times  nominal  damping.  As  shown  in  Figure  5-1,  the  optimum  damping  level 
produced  an  increase  in  speed  of  as  much  as  38  mph  on  one  of  the  smoothest 
courses  and  as  little  as  0.3  mph  on  one  of  the  roughest  courses.  This 
implies  that  the  nominal  damping  level  is  a  good  design  value  for  the 
rougher  terrain. 

Figure  5-2  shows  the  optimum  damping  factor  and  its  least  squares  regression 
line  as  a  function  of  terrain  RMS.  If  the  equation  of  the  regression  line 
were  used  to  predict  the  optimum  damping  factor  from  the  terrain  RMS,  then 
the  maximum  error  would  be  0.38,  the  average  would  be  0.12,  and  the  root 
mean  square  (RMS)  error  would  be  0.18.  This  amount  of  error  (especially  the 
0.38  maximum  value)  could  result  in  a  rougher  vehicle  ride  over  some 
terrains  than  using  the  nominal  damping  level. 

In  order  to  find  a  better  predictor  of  the  optimum  damping  level,  the 
terrains  were  frequency  analyzed  using  digital  filter  techniques  as 
described  in  Section  5.2.  of  this  report.  The  RMS  of  each  terrain  attribu¬ 
table  to  various  wavelengths  is  shown  in  Table  5-6.  This  type  of  terrain 
analysis  could  be  done  in  real-time  if  the  terrain  elevation  could  be  sensed 
and  fed  into  an  on-board  microcomputer.  Since  the  LK10934  concept  was 
designed  to  have  a  bounce  natural  frequency  of  68  cycles  per  minute,  there 
will  be  a  diminishing  effect  of  those  wavelengths  which  require  more  than 
0.88  second  for  the  vehicle  to  traverse  to  a  given  speed.  Thus,  an 
"effective  RMS*'  value  was  computed  for  each  terrain  as  the  square  root  of 
the  sum  of  the  squares  of  all  the  RMS  values  in  Table  5-6,  summing  from  zero 
to  the  first  frequency  band  whose  wavelength  mid-point  in  feet  exceeded  1.0 
times  the  vehicle  speed  in  ft/sec.  This  square  root  of  the  sum  of  squares 
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technique  is  used  instead  of  simple  summing  because  of  the  way  in  which  the 
digital  filtering  distributes  the  amplitude  of  a  signal  whose  frequency  s 
near  the  edge  of  one  of  the  wavelength  bands.  Figure  5-3  shows  the  optimum 
damping  factor  and  its  least  squares  regression  line  as  a  function  of  the 
speed  normalized  or  "effective"  RMS.  This  regression  line  can  predict  the 
optimum  damping  factor  with  a  maximum  error  of  0.19,  average  error  of  0.10, 
and  RMS  error  of  0.13.  This  improved  level  of  predictability  compared  to 
t,rraln  RM  that  ther.  la  a  d.flnlte  potential  for 

determining  the  optimum  damping  level  as  a  function  of  terrain  roughness  an 
frequency  content. 

The  lengths  of  the  nine  terrains  used  for  these  investigations  were  from  300 
to  500  ft.  If  the  roughness  level  or  frequency  content  of  a  terrain  changes 
significantly,  then  the  optimum  damping  level  also  changes.  A  measure  of 
the  variability  of  the  nine  aa:mple  terrains  was  obtained  by  computing  the 
elevation  RMS  over  several  200-foot  long  segments  of  a  terrain  and  then 
dividing  the  maximum  RMS  value  by  the  minimum  RMS  value  obtained.  The  indi¬ 
vidual  RMS  values  were  calculated  from  1  to  200  feet,  11  to  210  fee  ,  - 
220  feet,  and  so  on,  until  the  end  of  the  course.  If  the  course  is  uni¬ 
formly  rough,  then  the  variability  factor  computed  by  this  method  would  be 
l!o.  For  th^  nine  terrains  used  here,  the  variability  factor  shown  in  Table 
5-6  ranged  from  1.10  to  2.29  with  an  average  of  1.50.  Thus,  it  seems  that 
continuously  adjusting  the  damping  level  as  a  function  of 
roughness  and  frequency  while  traversing  a  course  should  all®”  Ims  L 
speeds  than  using  any  single  damping  level.  However,  to 
the  computer  would  require  considerable  effort  including  ^ 

digital  filtering  techniques  and  other  features  into  the  VRIDK  comp 
program.  This  has  not  yet  been  done. 

5.5.  Damping  Control  Baaed  on  Pitch  Rate 

Because  of  the  problems  involved  with  on-board  terrain  measurement  and  its 
prSent  eLly%evelop«ent  state,  It  was  decides  to  Investigate  the  potential 
for  ride  Improvement  based  only  on  a  knowledge  of  the  vehicle  motions.  In 
particular,  the  pitching  motions  of  the  vehicle  are  readily  measureable  with 
Ldern  instrumentation.  In  1975,  Dr.  R.  Beck  and  R.  Salemka  (Reference  5l 
reported  success  in  simulating  the  reduction  of  average  pitch  rate  by 
rapidly  switching  damping  blow-off  level  between  1.5 
nal  as  a  function  of  pitch  direction  (that  is,  up  or  down).  , 

Sfec?  on  Zorbed  poSar  was  not  rsportad.  It  Is  shown  below  that  e^derate 

success  In  reducing  absorbed  power  can  be  achieved  '>1  * 

slow  but  continuous  change  in  damping  as  a  function  of  the  time  his  y 
vehicle  pitch  rate. 

Computer  runs  with  some  additional  output  variables  implemented  VRIDE 
have  shown  that  instantaneous  pitch  rate  fluctuates  "«®’’  ”°^®. 
is  desirable  for  a  system  control  variable  or  for  ^^® 

of  the  terrain  (See  Figure  5-4).  The  following  formula,  which 
a  first  order  lag,  was  added  to  VRIDE  to  produce  a  running  average  value  of 

pitch  rate:  ,  ,,  i  • 

e_  =  lej  +  (T^N  -  1)  0a 
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where  9  is  the  instantaneous  pitch  rate,  '^a  is  the  running  average 
(filtered)  pitch  rate,  Tg.is  the  desired  time  constant,  and  N  is  the  number 
of  times  per  second  that  9  is  monitored  and  6  is  updated.  After  looking  at 
0Q  resulting  from  various  values  of  in  VRIDE  runs,  it  was  decided  to  use 
Tq  =  seconds.  This  caused  the  amplitude  of  the  pitch  rate  9^ 
by  almost  a  factor  of  2  compared  to  the  unfiltered  signal.  (Figure  5-5). 

It  was  decided  to  adjust  the  damping  rate  adjust  orifice  area  of  the 

damper)  as  a  function  of  averaged  pitch  rate,  85,  in  the  general  manner  as 
sketched  below. 


A  relatively  large  number  of  Iterations  were  required  to  obtain  an  adaptive 
suspension  control  scheme  of  the  type  outlined  above  which  would  improve  the 
ride  of  the  LKIOgS^*  on  moat  terrains  and  not  appreciably  worsen  it  on  any. 
Optimizing  this  control  scheme  requires  determining  the  optimum  values  of 
four  variables:  the  minimum  and  the  maximum  damping  levels,  and  the  slope 
and  the  x  axis  intercept  of  the  straight  line  which  transitions  from  the 
minimum  to  the  maximum  value.  After  many  iterations,  the  final  values  cho¬ 
sen  for  these  parameters  are:,  minimum  value  .667,  slope  1.667,  y  intercept 
0.0,  and  maximum  value  4.0.  9^  was  measured  in  rad/sec.  The  results 

obtained  for  the  LK10934  concept  vehicle  are: 

6  Watt  Speed  (MPH) 

Terrain  RMS  Standard  Suspension  Adaptive  Suspension 


APG-37  .68  51 .0 
APG-9  1.04  27.4 
APG-11  1.32  28.2 
APG-29  2.17  19.8 
APG-31  3.88  13.1 


62.3 

43.0 

33.2 
20.0 

14.3 


This  control  scheme  was  rudimentary  in  that  it  worked  with  only  the  single 
input  variable  9a.  Pitch  acceleration  or  road  arm  position  might  also  be 
very  useful  control  parameters.  Nevertheless,  substantial  increases  were 
demonstrated  in  the  vehicle’s  6-watt  speed.  Better  results  are  likely  with 
more  sophisticated  control  laws. 


The  relatively  small  improvements  produced  by  the  current  control  scheme  on 
the  higher  RMS  terrains  seem  to  be  due  to  a  requirement  for  high  damping  on 
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these  terrains  even  when  the  pitch  rate  is  temporarily  low.  Once  lowered  by 
the  control  scheme,  damping  is  not  increased  until  after  pitch  rate  has 
increased  and  thereby  already  adversely  affected  ride  quality.  Variable 
damping  as  a  function  of  spring  deflection  is  already  available  for  automo¬ 
bile  shock  absorbers.  A  few  VRIDE  runs  were  made  trying  out  this  technique, 
but  no  immediate  success  was  achieved.  However,  damping  variation  as  a 
function  of  spring  deflection  in  combination  with  pitch  rate  seems  like  a 
promising  area  for  future  investigation. 


The  amount  of  effort  required  to  obtain  the  adaptive  control  scheme  for  the 
LKIO934  points  out  the  need  for  a  more  organized  approach  or  a  well-defined 
method  for  adaptive  control  law  design.  For  some  particular  terrain  with  a 
relatively  constant  RMS  value  over  its  entire  length,  the  6-watt  speed  of 
the  vehicle  could  be  determined  using  VRIDE  with  different  constant  damping 
levels.  Enough  runs  would  be  made  to  determine  the  optimum  value  for  fixed 
damping.  This  process  could  be  repeated  for  terrains  covering  the  0.5  to 
4  0  inch  RMS  range.  The  average  value  of  the  running  average  pitch  rate 
would  be  recorded  for  each  run.  When  the  entire  process  was  completed, 
there  would  be  available  a  table  of  optimum  constant  damping  factor  versus 
terrain  RMS  and  average  pitch  rate.  From  this  table,  the  parameters  could 
be  determined  for  the  type  of  control  scheme  applied  to  the  LK10934  simula¬ 
tion.  Since  a  terrain  RMS  is  typically  not  constant  over  its  entire  length, 
and  since  the  four  second  running  average  pitch  rate  has  been  observed  to 
fluctuate  from  0.3  to  2.0  times  its  average  value  during  the  course  of  a 
single  VRIDE  run,  it  is  expected  that  fine  tuning  of  the  control  scheme 
would  still  be  required. 
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Figure  5-3.  Optimum . Damping  Factor  Vs.  Speed  Normalized  RMS 


5.6.  Application  of  On-Off  Passive  Dampers  To  a  Multiple-Wheeled  Vehlole 

This  section  presents  the  simulation  results  obtained  from  the  application 
of  the  concept  of  on-off  passive  dampers  (Reference  6,7)  to  a  ten-wheeled 
test  vehicle.  It  is  shown  that  improvement  in  the  vehicle  ride  performance 
over  test  terrain  is  possible  using  the  on-off  passive  dampers  in 

lieu  of  the  conventional  dampers. 

For  continuity  in  the  presentation,  a  simple  vibration 
a  conventional  damper  is  described  below.  The  on-off  passive  damper  is 

introduced  afterward. 

The  concept  of  the  on-off  passive  damper  is  then  applied  to  the  centered 
cLter  o/gravity  (eg)  configuration  of  the  ten-wheeled  HIMAG  test  bed 
veMoL!  ?he  perforLno.  of  the  on-off  peesLe  hamper  vibration  tnolatlon 
sJleine  was  simulated  in  the  VRIDE  computer  model.  The  results  are  doou- 
merited  below. 

The  differential  equations  of  motion  of  a  conventional  suspension  system 
Sat  is  fS  a  mass,  m  (acted  on  by  a  linear  spring  and  damper)  is  given  by: 

mx  =  -  b  (x-Xq)  -  k(x-xo) 

where  k  and  b  are  respectively  the  spring  constant  and  the  damping  coefficient. 

_ j _ A.  •  ^  i  \  ■!  «  vrAKi  Kv  • 


(2) 


The  values  of  k  and  b  are  usually  chosen  to  provide  the  desired  damping 
ratio  or  the  desired  resonant  frequency  ^n. 

An  operational  diagram  and  a  control  ]ogic  for  a  simple 

suspLsion  system  (References  1  &  2)  are  / 

oontrol  variable  of  the  system  is  the  adjustable  damper  ^ 

Sr llble  damper  is  assumed  to  be  controlled  by  a  high  band-width  damper 

orifice  actuator. 

The  equation  of  motion  for  the  on-off  passive  damper  system  shown  in  Figure 
5-6  is  given  by: 

mZ  =  k  (7-Zo)  -  b(t)  (Z-Zq). 

As  shown  in  Figures  5-6  and  5-7,  the  coefficient  b(t)  takes  on  the  constant 
talue  b  or  0  depending  on  whether  Z  and  (Z-Zq)  are  of  the  same  sign  or  not. 

That  is: 


0m 

n 

1 

- 

• 

X 

0  1 

0  0 

Xo 

.  1 

+ 

I 

• 

-k/m  -b/m_ 

■fk/m  •♦b/mj 

Xo_ 

4  • 


b(t)  = 


b,  for  Z(Z-Zo)>0; 
0,  for2(Z-7o)<0- 


(4) 
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A  state  space  equation  for  (3)  nay  be  written  as: 

z  r  0  1  z  0  0 

zj  (-k/m  -b/mj  [  zj  [k/m  b(t)/m^ 

The  on-off  passive  suspension  system  defined  by  equations  (3)-(5)  has  been 
shown  in  References  6  and  7  to  provide  a  desirable  vibration  isolation 
characteristic  which  filters  out  high  frequency  excitation  inputs  when 
applied  to  a  mass  supported  by  a  single  spring  and  damper.  In  this  report, 
the  logic  described  by  eq.  (4)  is  applied  to  a  ten-wheeled  vehicle  which  is 
simulated  using  a  modified  version  of  the  VRIDE  program  (Reference  3)  on  the 
PRIME  computer  facility  at  the  US  Army  Tank-Automotive  Command  in  Warren, 
Michigan.  In  the  simulation,  the  logic  of  eq.  (4)  is  installed  indepen¬ 
dently  for  each  individual  wheel  suspension  subsystem.  The  actual  parame¬ 
ters  of  the  center  eg  configuration  of  the  HIMAG  test  bed  vehicle  were  used. 
The  VRIDE  model  has  been  extensively  validated  against  field  test  data  for 
this  vehicle. 

For  purposes  of  suspension  system  evaluation,  this  vehicle  was  simulated 
running  over  a  set  of  five  Ft.  Knox  test  courses  which  represent  a  range  of 
quite  smooth  to  very  rough  cross-country  conditions.  The  vertical  eleva¬ 
tions  of  these  courses  were  measured  at  1-foot  horizontal  intervals.  These 
values,  with  linear  interpolations  between  data  points,  were  used  as  a 
forcing  funct^.on  at  the  bottom  of  the  roadwheels  in  the  VRIDE  model. 


For  comparison  purposes,  the  simulation  runs  were  carried  out  with. 

(I)  the  vehicle  having  conventional  dampers; 

(II)  the  vehicle  having  the  on-off  passive  dampers  (logic  of  equation  4); 

(III)  the  vehicle  having  the  on-off  passive  dampers  (logic  of  equation  4) 
but'  with  increased  damper  coefficient  of  1.5  tiroes  its  original 
values . 

In  the  simulation,  these  vehicles  were  "test  driven"  over  the  simulated 
terrains  at  various  constant  speeds.  The  performance  of  a  vehicle  is  moni¬ 
tored  at  speeds  of  5,  10 . 70  mph..  The  list  p.f  monitored  variables 

includes  the  absolute  values  of  0  max,  0  max,  and  0  max  (maximum  pitch 
angle,  velocity,  and  acceleration);  0rms,  0rms,  and  0rms  (RMS  pitch  angle, 
velocity  and  acceleration),  yrms  and  :^ax  (RMS  and  maximum  vertical  velocity 
at  eg);  jimax  (maximum  vertical  acceleration  at  eg;  y  max  and  y  RMS 
(maximum  and  RMS  vertical  acceleration  of  the  driver)  and  the  driver- 
vertical-absorbed  power,  Pd. 

A  series  of  VRIDE  runs  was  made  comparing  the  on-off  damper  (II)  to  the 
standard  HIMAG  damper  (I)  on  the  five  Fort  Knox  terrains.  Improvement  in 
driver-absorbed  power  with  the  on-off  dampers  was  shown  on  all  the  courses. 

On  the  average,  for  any  given  speed,  the  driver-absorbed  power  was  reduced  by 
almost  50)1  on  the  smoother  courses  and  by  at  least  on  the  rougher  cour¬ 
ses.  The  speed  at  which  6  watts  was  reached  increased  by  more  than  20  mph 

on  the  smooth  courses  and  less  than  1  mph  on  the  roughest  courses. 

For  the  median  roughness  course  (.173  in.  RMS),  plots  were  made  of  several 
vehicle  outputs  vs.  speed.  These  show  that  the  two  damping  systems  per¬ 
formed  basically  the  same  in  terms  of  maximum  pitch  angle,  maximum  pitch 
rate,  and  RMS  pitch  angle.  The  on-off  damper  was  somewhat  better  in  terms 
of  maximum  pitch  acceleration,  RMS  pitch  rate,  and  maximum  vertical  acce¬ 
leration  of  the  driver,  and  in  RMS  vertical  velocity  of  the  eg.  In  terms 
of  maximum  and  RMS  vertical  acceleration  of  the  eg  and  RMS  pitch  accelera¬ 
tion,  the  on-off  damper  was  definitely  better,  at  times  by  as  much  as  30%. 
When  the  on-off  damper  was  better,  it  was  consistently  better  for  speeds  in 

the  5  to  70  mph  range,  (up  to  12  watts).  The  driver-absorbed  power  vs. 

speed  plots  are  given  below: 


1 

I 


<?lnce  the  on-off  damper  system  (II)  was  just  added  on  to  an  already  designed 
rtamDinir  svstem  (I),  it  seems  that  even  better  performance  might  be  obtained 
K  the  on-off  effect  was  considered  in  the  overall  suspension  design  pro¬ 
cess  Since  the  on-off  damper  is  "on”  only  when  it  tends  to  Improve  the 
ride’  it  seems  that  a  higher  level  of  damping  might  improve  the  ride  even 
more’  This  might  also  give  relatively  more  improvement  on  the  rougher  cour¬ 
ses  where  it  has  not  yet  been  achieved.  For  these  reasons,  damping  system 
(III)  was  formulated  with  the  damping  level  1.5  times  the  nominal  level  an 
with  the  on-off  control  law.  The  HIMAG  vehicle  was  simulated  over  test 
course  FTKN910  (2.2M  in.  RMS). 

For  the  test  runs  over  this  simulated  terrain  with  the  standard  damping,  the 
vehicle  exhibits  a  distinct  resonance  characteristic  over  the 
25  mph  for  driver-absorbed  power;  maximum  pitch  angle,  pitch  rate,  and  pitch 
acceleration;  and  RMS  of  pitch  angle,  pitch  rate,  and  pitch  acceleration. 

Also  maximum  vertical  acceleration  at  the  driver  and  eg  locations  ® 

very ’strong  peak  in  the  15-25  mph  range.  The  on-off  damper  system  with  the 
1.5^nominal  damping  level  shows  a  marked  decrease  in  the  resonance  peak 

all  these  variables. 

In  aeneral  with  on-off  dampers,  improvements  in  the  simulated  performance 
ff  ?he  v^ilc^e  are  obtalned%n  all  the  monitored  v-iables  when  compared  to 
standard  dampers..  Significant  .reduction  In  .the 

are  obtained  ln|fe|max,  0  RMS,  YRMS,  Y  RMS,  |Y|max,|  Y^^j  maxjv^l  RMS,  and  Pf,. 
Superior  high  speed  performance  area  achieved  in  terms  of  considerably 
reduced  6'  RMS,  Y  RMS,  Y  RMS,  and  Pd  (driver-absorbed  power). 

For  the  simulated  vehicle  with  on-off  passive  dampers  having  Increased 
damper  cS?icients  of  1.5  times  its  original  values  excellent  reduction  xn 
thfresonance  transmittivity  is  found  in  all  the 

coLared  to  the  two  previous  cases.  The  high  speed  characteristics  for  the 
increased  damping  system  are  similar  to  those  of  on-off  passive  dampers 
kse  II.  Due  to  the  increased  damping,  very  slight  f 

speed  as  well  as  low  speed  performance  may  be  noticed  in  |Y  RMSj  , |Y  max/ 

|Y§max  I ,  Iy^rms ( ,  and  Pd.  However,  the  decrease  in  the  resonance  peaks  in  t 
1 5-25  mph  range  are  much  more  noteworthy . 

For  clarity  of  presentation  the  improvement  in  the  ride  performance  obtained 
from  the  sLula?ion  results  are  summarized  in  Table  5-7.  Approximate  improve¬ 
ment  indices  are  included  for  ease  of  comparison. 

The  work  described  in  this  section  is  ongoing  at  Oakland  University  at  the 
time  of  this  report.  However,  it  appears  quite  promising.  The  complete 

at  4  greater  level  of  detail  In  a  fortheonlng 

contractor  report. 
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Table  5-7 .  Improvement  Chart  of  Ride  Characteristics  for  Vehicle  with  On-Off  Passive  Dampers 
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moderate  improvement  (approximately  5-25%) 
significant  improvement  (approximately  25-50%) 
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APPENDIX  A 


PRELIMINARY  STUDY  FOR  ADAPTIVE 
SUSPENSION  SYSTEM  PROGRAM 


A-1 


Freliininary  Study  for 
Adaptive  Suspension  System  Program 


Prepared  by:  Francis  B.  Hoogterp 
Applied  Rsch  Func 
Tk-Autmv  Concepts  Lab 
September  1980 


An  overall  Adaptive  Suspension  System  Program  Summary  is 
included  as  Appendix  I.  Certain  general  requirements  of  this 
system  are  of  particular  concern  here.  Since  the  system  must 
be  capable  of  operating  at  speeds  up  to  50  mph  and  take  pro¬ 
file  measurements  at  least  every  foot,  this  imposes  a  minimum 
rate  requirement  of  74  samples  per  second.  It  should  be 
noted  also  that  at  10  mph  the  vehicle  traverses  only  73  feet 
in  the  5  seconds  allotted  for  system  start  up. 

The  complete  sensor-processor  design  cannot  be  postulated 
until  the  results  of  the  suspension  optimization  effort  have 
been  analyzed  and  the  decision  is  made  as  to  the  manner  in 
which  the  suspension  system  will  be  alterable.  Much  can  and 
has  been  accomplished  in  the  area  of  sensor  specifications 
and  processing  requirements. 

Since  we  hope  to  use  the  results  of  the  T95  suspension 
optimization  study  to  formulate  the  basis  for  our  control 
choices,  we  have  automatically  constrained  our  system  to 
process  terrain  profile  data  to  obtain  rms  information  and/or 
the  frequency  content, of  the  terrain.  The  traditional  method 
for  characterizing  terrain  roughness  is  to  calculate  the  rms 
(root  mean  square;  of  the  profile.  This  is  done,  however, 
only  after  passing  the  profile  measurements  through  a  high 
pass  filter  with  a  cutoff  frequency  of  1/60  cycle  per  foot. 
This  characterization  ignores  differences  in  the  slopes  of 
the  terrain  PSD's  which  depict  the  varying  predominant 
frequencies  which  are  truely  present  in  different  profiles. 

One  approach  to  processing  the  profile  data  would  be  to 
use  an  FFT  (Fast  Fourier  Transform)  algorithm  to  estimate  the 
PSD.  This  approach  has  several  drawbacks.  An  FFT  on  300  or 
400  points  can  take  considerable  computation  time.  Also  each 
time  it  is  to  be  updated,  the  computation  must  be  restarted 
and  all  the  data  points  must  be  available  simultaneously. 

Another  approach  would  be  to  take  the  digitized  signal 
and  run  it  through  several  high  pass  and  low  pass  filters 
simultaneously  (see  Figure  1). 


Figure  1 
A-3 


The  output  of  each  filter  would  be  squared  and  averapd 
to  Rive  the  mean  square  value  of  the  filtered  signal.  The 
mean  square  outputs  can  then  be  differenced  to  obtain  the  mean 
square  contribution  between  each  two  adjacent  cutoff  frequencies. 
The  averaging  of  the  mean  squared  signal  can  be  performed  with 
a  time  constant  selected  to  average  the  most  recent  4  or  5 
seconds  of  data.  In  this  way  the  data  need  not  be  stored  and 
the  output  is  continuously  available.  In  other  words  the 
averaging  need  never  be  restarted  in  order  to  update  the 
information  required  to  control  the  suspension  system. 


This  latter  approach  is  definitely  recommended  for  imple¬ 
mentation.  The  number  of  filters  required  will  probably  be 
four  unless  an  additional  one  is  required  due  to  high  frequency 
noise.  The  recommended  configuration  is  shown  in  Figure  L. 


A  couple  of  final  notes  should  be  made.  Since  the  effort 
is  to  define  the  profile  independently  of  the  vehicle  velocity 
the  accelerometer  signal  should  be  digitized  at  a  constant 
step  size  in  feet.  This  can  be  accomplished  by  instrumenting 
the  odometer  to  output  a  signal  perhaps  every  1/16  ft.  of 
horizontal  travel.  This  signal  can  then  be  used  directly  or 
counted  down  to  obtain  the  digitizing  signal  remembering  of 
course  that  the  maximum  sample  distance  acceptable  is  1  toot. 
Also  care  must  be  taken  to  insure  that  the  processor  can  keep 
up  in  the  worst  case  situation  (i.e.  at  50  mph  and  with  control 
changes  to  be  made) . 

The  filters  described  can  best  be  implemented  digitally. 

The  transformation  to  the  spatial  domain  is  then  achieved  simply 
by  maintaining  a  constant  sample  interval  in  feet.  If  the  fil¬ 
ter  were  implemented  with  analog  components,  either  integration 
gains  or  cutoff  frequencies  would  have  to  vary  with  the  vehicle 
velocity.  This  would  greatly  complicate  the  design. 


The  adaptive  suspension  system  implementation  will  defin¬ 
itely  require  a  16  bit  microprocessor.  The  implementation  of 
the  filters  should  be  done  with  integer  arithmetic  and  will 
definitely  require  an  integer  arithmetic  logic  unit  . 

The  program  will  be  stored  in  EPROMS  and  128  bytes  of  RAM  for 
working  memory  should  suffice.  If  the  acceleration  is  sampled 
each  1/4  foot  then  at  50  mph  the  maximum  time  available  could 
be  3.4  msec  to  complete  all  computations  in  the  loop.  With 
this  in  mind,  the  filters  should  probably  be  only  second  order. 
Also  if  it  becomes  necessary,  high  speed  16  bit  by  16  bit  multi 
pli  ers  are  available. 


A  diagram  of  the  probable  microprocessor  layout  is  shown  in 
Figure  3.  The  time  constant  for  the  averaging  circuit  is  1/T 
and  should  probably  be  selected  to  be  about  3  seconds.  Hope¬ 
fully  the  first  order  averaging  circuit  will  be  adequate. 
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FIGURE  3  Microprocessor  Configuration 
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Sensor  Alternatives 


Three  specific  types  of  sensors  were  considered  for 
incorporation  into  the  adaptive  suspension  system.  These  were 
(1)  a  sonar  reflective  sensor,  (2)  a  pressure  transducer,  and 
(3)  an  accelerometer.  Each  of  these  sensors  has  certain  ad¬ 
vantages  and  disadvantages  which  will  be  discussed  briefly. 

Sonar  Reflective  Sensor 

A  sensor  of  this  type  could  be  mounted  to  the  underside  of 
the  hull,  perhaps  near  the  center,  and  could  provide  the  dis¬ 
tance  between  the  hull  and  the  ground.  The  particular  unit 
looked  at  (made  by  Polaroid)  was  capable  of  measuring  distances 
of  from  0.9  ft.  upt  to  35  ft.  to  the  nearest  0.1  feet.  The 
sampling  frequency  was  limited  to  5  readings  per  second.  This 
maximum  sampling  frequency  was  designed  into  the  sensor  in  , 
order  to  obtain  a  maximum  range  of  35  feet.  The  sound  requires. 
62  ms  to  travel  the  35  feet  out  and  back  and  therefore,,  limits 
the  design  frequency  to  something  less  than  16  readings  per - 
second,  far  less  than  required  in  this  application.  It  the 
range  were  limited  to  say  3.5  ft  maximum,  the  the  maximum 
transmission  time  would  be  reduced  to  6.2  ms  and  perhaps  a 
rate  of  100  samples  per  second  could  be  achieved.  The  original 
requirement  was  to  sample  at  least  73  times  per. second  (every 
foot  at  50  mph).  It  should  be  noted  that, if  a  radar  device  were 
used  a  50  nanosecond  output  pulse  would  cause  a  minimum  range 
of  25  ft.  and  therefore  would  not  be  practical  for  this 
application. 

Another,  and  perhaps  more  serious  problem  exists  with  the 
use  of  any  reflective  type  sensor.  The  information  received 
would  provide  only  the  relative  distance  between  the  hull  and 
ground  at  any  point  in  time.  This  is  not  the  same  as  providing 
the  terrain  profile  height.  A  detailed  analysis  was  performed 
to  investigate  the  degree  of  correlation  between  this  measured 
signal  and  the  actual  profile  height.  The  T95  was  simulated 
over  three  terrains  of  various  surface  roughness  (i*15  in,  l.\l 
in,  and  3.88  in  rms)  and  at  various  vehicle  velocities. 

The  resulting  correlation  between  the  ground  profile  and 
the  T(ftletive  displacement  between  the  hull  and  the  terrain  can 
be  seen  in  the  following  Figures  4  thru  15. 

Each  of  the  profiles  wa’s  filtered  with  a  two  pole  high  pass 
filter  with  a  spatial  cutoff  frequency  of  1/60  cycle/foot.  The 
two  signals  of  interest  diverge  greatly  at  times  indicating  the 
vertical  movement  of  the  hull.  Even  though  this  hull  movement 
has  a  lower  frequency  content  than  much  of  the  profile  signal, 
it  is  still  significantly  within  the  band  of  interest. 
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In  order  to  successfully  utilize  a  reflective  sensor  it 
would  also  be  necessary  to  include  instrumentation  to  determine 
the  hull  vertical  position  relative  to  the  earth  and  to  remove 
its  effect.  This  could  be  done  with  a  hull  mounted  accelerometer 
and  would  permit  the  estimation  of  the  roadwheel  height,  and 
therefore  the  profile  height  from  the  combination  of  the  two 
sensor  outputs. 

Pressure  Transducer 

It  was  conjecured  that  a  pressure  transducer  installed  within 
the  hydropneumatic  suspension  system  could  provide  a  reading 
that  could  readily  be  converted  to  spring  force.  The  processor 
could  then  make  the  conversion  through  a  table  look  up  to 
obtain  the  spring  deflection.  This  technique  would  provide  a 
relatively  clean  signal  and  the  sensor  would  not  be  subjected 
to  the  environment.  Also  the  samples  could  be  digitized  under 
command  of  the  microprocessor  and  would  not  be  restricted  to  the 
relatively  low  sample  rates  of  the  sonar  reflective  sensor. 
Unfortunately  it  would  also  only  provide  relative  information 
of  profile  height  with  respect  to  the  hull's  vertical  position. 
Therefore,  a  hull  accelerometer  would  also  be  required  in  order 
to  obtain  a  good  estimate  of  the  profile  height. 

Accelerometer 


The  third  sensor  option  considered  was  to  use  a  single 
accelerometer  mounted  on  the  third  roadarm  as  near  to  the  road- 
wheel  as  possible.  The  accelerometer  should  be  mounted  in  a 
position  such  that  a  true  vertical  reading  is  obtained  when 
approximately  2  inches  of  jounce  travel  have  been  encountered. 

The  2"  value  is  somewhat  arbitrary  but  is  selected  in  an  attempt 
to  minimize  the  geometric  errors  caused  by  the  roadarm  rotation. 

This  option  was  selected  despite  the  problems  of  mounting 
the  accelerometer  in  an  exposed  location  and  despite  the  added 
geometric  considerations  due  to  roadarm  rotation.  Since  an 
accelerometer  would  also  be  required  in  other  sensor  options 
(though  it  could  be  mounted  internally),  the  single  accelerometer 
approach  seems  more  reasonable  at  least  in  the  test  bed  phase. 
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Terrain 


APG30 


APG29 

APG31 


Roughness 

(rms) 

Simu lation 

MPH 

Results 

WATTS 

(1.15") 

30 

0.88 

35 

1.33 

40 

2.00 

45 

3.01 

50 

4.40 

(2.17') 

20 

3.97 

25 

5.89 

(3.88") 

12 

4.55 

15 

11.58 

Table  1 

T95  VRIDE  Simulation  Results 


ADAPTIVE  SUSPENSION  SYSTEM  PROGRAM  SUMMARY 

OBJECTIVE:  To  develop  a  suspension  system  with  the  inherent ^ 

capability  to  adjust  operating  characteristics  to  suit  changing 
terrain  inputs.  This  capability  would  result  in  improved 
vibration  isolation  with  accompanying  gains  in  ride  quality 
and  gun  platform  stability. 

GENERAL  OPERATION:  On-vehicle  sensors  will  measure  the  terrain 
profile  as  it  is  being  traversed.  The  suspension  characteristics 
needed  for  optimum  performance  in  this  terrain  will  be  determined 
from  a  suspension  optimization  program  in  the  on-borad  micro- ^ 
processor  memory.  The  suspension  system  will  then  be  automatically 
adjusted  to  provide  optimum  performance. 

HARDWARE  REQUIREMENTS : 

a.  Terrain  Sensors 

b.  On-Vehicle  Microprocessor 

c.  Servo  Controls  for  Suspension 

d.  Variable  Suspension  System 
SOFTWARE  REQUIREMENTS: 

a.  Microprocessor  Logic  Package 

(1)  Process  Terrain  Data 

(2)  Identify  the  terrain  range  in  which  the  vehicle  is 
operating. 

(3)  Determine  suspension  requirements  for  existing 
terrain  range,  (from  pre-established  program  in  microprocessor 
memory) . 

(4)  Send  Servo  control  signals  to  make  necessary  sus¬ 
pension  adjustments. 

b.  Suspension  Optimization  Hybrid  Computer  Program 

FUNCTIONAL  OPERATION:  Sensors  will  measure  the  terrain  profile 
a"s  Tt  ITs  being  traversed  and  provide  input  to  the  on-board  micro- 
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processor.  The  microprocessor  vill  process  the  terrain  sensor 
data  to  provide  a  measure  of  the  terrain  profile  such  as  the 
PSD  and/or  RMS.  This  measured  terrain  profile  data  will  then 
be  compared  with  a  set  of  representative  terrain  ranges  from 
the  hybrid  suspension  optimization  program  in  the  microprocessor 
memory  and  the  terrain  range  in  which  the  vehicle  is  operating 
will  be  identified.  The  corresponding  optimimum  suspension 
parameters  for  this  terrain  will  also  be  simulataneously 
identified.  The  microprocessor  will  then  send  the  required 
Servo  control  signal  to  make  the  necessary  suspension  adjustments. 

GENERAL  REQUIREMENTS: 

a.  System  must  be  capable  of  functioning  at  vehicle  speeds 
of  up  to  50  MPH. 

b.  System  must  be  capable  of  functioning  on  terrain 
profiles  up  to  A  RMS. 

c.  Initial  terrain  profiles  must  be  developed  in  not 
more  than  five  seconds  and  updated  every  second. 

d.  Terrain  profile  to  be  measured  in  one  foot  increments. 

e*  System  will  incorporate  a  manual  overdrive. 

PROGRAM  DEVELOPMENT  ELEMENTS: 

a.  Hybrid  computer  suspension  optimization  program 
development. 

b.  Terrain  sensor/microprocessor/microprocessor  logic 
package  development . 

c.  T95  vehicle  modification  and  evaluation. 
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FIGURE 


T95  Al  TON  TANK 

30  MPH  ON  APG30  -Cl. 15  IN.  RMS)  10  NOU  76 
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T95  -41  TON  TANK 

UHEEL  3  TO  CHASSIS  DISPLACEHENT 

HP  .7  HZ  UHEEL  NUMBER  3  PROFILE  ELEVATION 
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FIGURE 


T95  41  TON  TANK 

45  NPH  ON  APG30  (1.15  IN.  RNS )  10  NOV  76 
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FIGURE 


T95  41  TON  TANK 

WHEEL  3  TO  CHASSIS  DISPLACEMENT 

HP  1.1  HZ  WHEEL  NUMBER  3  PROFILE  ELEVATION 
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FIGURE  10 


T95  41  TON  TANK 
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figure  12 


T95  41  TON  TANK 

12  flPH  ON  APG31  (3.88  IN.  RNS )  16  NOU  76 
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